Weak correlation effects in the Ising model on triangular-tiled hyperbolic lattices 
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The Ising model is studied on a series of hyperbolic two-dimensional lattices which are formed 
by tessellation of triangles on negatively curved surfaces. In order to treat the hyperbolic lattices, 
we propose a generalization of the corner transfer matrix renormalization group method using a 
recursive construction of asymmetric transfer matrices. Studying the phase transition, the mean-field 
universality is captured by means of a precise analysis of thermodynamic functions. The correlation 
functions and the density matrix spectra always decay exponentially even at the transition point, 
whereas power law behavior characterizes criticality on the Euclidean fiat geometry. We confirm 
the absence of a finite correlation length in the limit of infinite negative Gaussian curvature. 
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I. INTRODUCTION 

An increasing interest in the thermodynamic behav- 
ior of various physical models on non-Euclidean (curved) 
surfaces has been persisting for about two decades, due 
to recent experimental fabrication of soft materials with 
conical geometry [l[ and magnetic nanostructures which 
exhibit negatively curved geometries [HI]. Curved ge- 
ometries are also relevant in the theory of quantum grav- 
ity [B|, Hj]. In this context, several statistical models have 
been investigated on simple negatively curved geome- 
tries, such as the Ising model LZ H9L the g-state clock 
models [MED], and the XY-model [1J] . 

A typical example of the negatively curved geometry 
is represented by the two-dimensional discretized hyper- 
bolic surface (lattice) which is characterized by a con- 
stant negative Gaussian curvature. Among the varieties 
of lattice surfaces, we choose, for simplicity, a group of 
regular lattices that are constructed as tiling of congruent 
polygons of the p-th order with the coordination num- 
ber q. On the hyperbolic (p, q) lattices, the relation 
(p — 2)(q — 2) > 4 is satisfied, in contrast to the rela- 
tion (p — 2)(q — 2) = 4 on the Euclidean flat geometry. 
Figure [1] shows two examples, the (3,7) and (3,13) lat- 
tices where the whole lattice is mapped onto the Poincare 
disk [13]. 

In general, the number of the lattice sites within a 
certain area increases exponentially with its diameter on 
such hyperbolic lattices. This exponential increase limits 
efficiency of numerical studies of statistical models, such 
as the Ising model on the (p, q) lattice. In particular, 
applications of Monte Carlo simulation face difficulties 
in the scaling analysis around the phase transition. Also 
transfer matrix diagonalization can not easily be applied 
due to the non-triviality in the construction of the row- 
to-row transfer matrices. 

Despite these difficulties, one can evaluate the parti- 
tion function by means of Baxter's corner transfer ma- 
trix formalism [l4| even for the hyperbolic (p, q) lat- 
tices. In this article, we use a flexible numerical im- 




FIG. 1: Poincare disk representation of the hyperbolic lattices 
created by triangular tessellation p — 3 with the coordination 
numbers q = 7 (left) and q — 13 (right). 



plementation of Baxter's method, so-called the Corner 
Transfer Matrix Renormalization Group (CTMRG) al- 
gorithm, which has been used as a tool in the computa- 
tion of the partition function for (flat) t wo- and three- 
dimensional classical s pin systems [15l - fl7j . In our previ- 
ous reports [^, [TO, fl8l Il9| we considered the hyperbolic 
(p. q) lattices, typically for the case with q = 4, where 
the whole lattice can be divided into four quadrants, the 
'corners'. For the Ising model on the (p, 4) lattices, the 
mean- field universality was found 0, H| ■ 

The hyperbolic (p, q) lattice with an arbitrary coordi- 
nation number q other than four has not been addressed 
by use of the CTMRG method yet. For this case the nu- 
merical renormalization procedure of the corner transfer 
matrices requires a technical extension upon the estab- 
lished numerical procedure for the (p, 4) lattices. In this 
article we introduce a new procedure which is valid for 
general values of q and find the thermodynamic prop- 
erties of the Ising model on a wider class of the (p, q) 
lattices. In particular, the triangular tessellation (p = 3) 
and the coordination number q > 6 are investigated as 
representative examples. 

This article is organized as follows. In Sec. II we define 
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the Ising model on the (p, q) lattices. In Sec. Ill the recur- 
rent renormalization algorithm of the CTMRG method 
is introduced. The application of CTMRG to the (3, q) 
lattices is explained starting from q = 6 and increasing q. 
Numerical results on the spontaneous magnetization and 
energy are presented in Sec. IV, with a detailed analysis 
of the g-dependence of the phase transition temperature 
and the corresponding scaling exponents. In Sec. V the 
quantum entropy and the scaling behavior of the correla- 
tion functions are observed. We also analyze the effects 
of the Gaussian curvature on the correlation length. We 
summarize the result in the last section. 



II. THE LATTICE MODEL 



Consider the Ising model with the Hamiltonian 



H{a) = - J ^ (TjCTj - hyVj 

{id} {i} 



(1) 



defined on the hyperbolic (p, q) lattices. We here use the 
standard notation (p, q) where the first integer p corre- 
sponds to the regular polygons with p sides (or vertices) 
and where the second one q stands for the coordination 
number which is the number of polygons meeting in each 
vertex. Throughout this article we focus on the trian- 
gular tiling on the (3,q) lattices only. The Ising spin 
variables a i =t or J, are located on the vertices. The 
first term in H(a) represents the ferromagnetic coupling 
(J > 0) between the nearest-neighboring Ising spins a i 
and <jj , and the second represents the effect of the exter- 
nal magnetic field h. Then, the partition function 



Z = ^2 exp 

M 



k B T 



(2) 



is given by the sum of the Boltzmann weights over all 
spin configurations which are denoted by {<r}. Here, k B 
and T, respectively, are the Boltzmann constant and the 
temperature. 

On any (p, q) lattice, the Boltzmann weight of the 
whole system can be represented as the product of the 
local Boltzmann weights attributed to the particular 
p-gons. In this study we define the local Boltzmann 
weights which are consistent with the triangular tessel- 
lation (p = 3). For a reason we explain in the follow- 
ing, each local Boltzmann weight W B is constructed by 
a pair of adjacent triangles cr a a b <Jd and a b a c ad as shown 
in Fig. 2. The local Boltzmann weight Wb for this pair 
of the triangles is then given by 



W B (a a a b a c <7 d ) = exp 

/ 

+&d&a + ZObOd) 



J 



2k B T 



(a a a b + o b a c + o c Od 



qk B T 



(a a + 2er& + <T C + 2er d ) 



(3) 



The factor 2 of 2<j b a d arises from the fact that a b and a c 
are shared by two triangles, under the tessellation of "bi- 
triangular" Boltzmann weights. Also the factor 2 appears 



in 2<76 and 2ad since the effect of external magnetic field 
h should be counted for both upper and lower triangles. 
Under these factorizations, we proceed the calculation by 
the CTMRG method HEl]. 

The standard numerical formalism based on the di- 
agonalization of the row-to-row transfer matrix is not 
easily applied under hyperbolic geometries. It has been 
shown that the CTMRG method works as an alterna- 
tive when the (p, 4) lattice is considered under the 
condition p > 4. Recall that the (p, 4) lattice can be di- 
vided into four equivalent quadrants by two perpendicu- 
lar geodesies, and it is easily understood that each quad- 
rant corresponds to the corner transfer matrix 0, [l8| . 
Such division of the whole system is not generally admis- 
sible for the (3, q) lattices, which is under our interest, 
in particular, when q is odd; we tackle this case in the 
following. 



III. RECURRENT RG SCHEME 

Let us consider the generalization of the CTMRG 
method to the (3,g) lattice with q > 6. The Boltzmann 
weight of the whole (3,q) lattice can be represented by 
the product of the q identical corner transfer matrices C q 
surrounding the central spin a. In this picture we can 
express the partition function in the product form 



(4) 



?i,«2,...,?,i=i 



where ^ and which appear as the parameters of 

the corner transfer matrix (^(cr^^j+i), are the block spin 
variables corresponding to chains of the spins from the 
central spin a towards the system boundary. We have 
assumed the cyclic order around cr, and thus = £i 
is satisfied. Throughout this paper we use the counter- 
clockwise index ordering for the spin variables included in 
corner transfer matrices Cq(cr^j^j + i), starting from any 
one of the two-state variables from <j a to a e i- Also the 
renormalized spin variables £j are aligned in the same or- 
dering, as shown on the red triangles in Figs. [U [31 and 

El 

Let us explain the recursive construction of the corner 
transfer matrix C g (cr^^+i) with respect to the successive 
area expansion of the whole system Q. For a tutorial 
purpose, we start from the (3, 6) lattice where the system 
is on the flat surface, and treat the cases q > 6 afterward. 

In contrast to the original CTMRG formulation [l5| , it 
is important to introduce two different kinds of 'half-row 
transfer matrices' C q and lZ q which are used for the area 
expansion of the corner transfer matrix. On the (3, 6) 
lattice, the area expansions of the transfer matrices Cq 
and TZq are performed as 

C 6 {a d cr a cr b a c ^2) = WB(cr a cr fc cr c cr (i )£ 6 (o' ti cr c ^i$2), (5) 
^(o-cOdCTaCr^i^) = W B (aa<J b a c a d )TZ 6 (a c a b ^ 2 ) (6) 



FIG. 2: (Color online) Graphical representation of the exten- 
sion process of the left transfer matrix L 6 , the right transfer 
matrix IZq, and the corner transfer matrix C% on the (3,6) 
lattice which are defined by Eqs. (©-(JTJl. The filled symbols 
correspond to the variables which have to be summed up. The 
two-state and multi-state variables, respectively, are denoted 
by a and £. 




FIG. 3: (Color online) The expansion process of £7, IZ7, and 
C 7 on the (3, 7) lattice. 



where the position of each spin variable is graphically 
depicted in Fig. [2j Similarly, the corner transfer matrix 
Cq are expanded as 

C§(VdOa°b£,l£,i) = ^ WB{0a0b<rcad)£>G{?dOc&Z,i) 

xC 6 (^26)ft6(^66k). (7) 

The recursive expansion procedure in CTMRG can be 
initiated by setting C 6 (a a a b a c a d ) = TZ 6 (a d a a a b a c ) = 
yVB(a- a crba- c a d ) and C 6 (a a a b a d ) = w B(c a cr&cr c crd) 
where the multi-spin variables £ are identical with the 
Ising ones a at the beginning. In the following, we do 
not write spin variables explicitly for book keeping. 

We now generalize the above-mentioned expansion pro- 
cess for the (3,(7) lattices, when q>7, where the hyper- 
bolic surface geometry is realized. Drawing the lattice, 
such as shown in Fig. 1, and analyzing the inner structure 
of the corners, one can derive a set of recursive relations. 
The g-dependent corner transfer matrix 

C q = Y, W B£ q C q K q (8) 

is a slight modification of Eq. ([7]). The relation is graph- 
ically shown in Figs. |3] and 2] for the two representative 
cases. Similarly, for the 'half- row transfer matrices', we 



FIG. 4: (Color online) The expansion process of £ 13 , 72. 13 , 
and C 13 for the (3, 13) geometry. 

obtain 

r, = ^w b c;"k,9 (10) 

where n q is the multiplicity of C q given by 

n q = ^—^ — = max |n G Z | n < ^—^ — j . (11) 

In contrast to Eqs. 5 and 6, the corner transfer matrices 
appear in the expansion process of C q and lZ q when q > 7. 
The extended transfer matrices C q , lZ q , and C q re-enter 
the right hand sides of Eqs. ([5])- (fT0|) . 

The expansion process successively increases the sys- 
tem size by expanding the matrix dimensions of C q , TZ q , 

and C q . To prevent the exponential grow of computa- 
tional effort, we introduce the density- matrix renormal- 
ization scheme 0, UH. Let us express the block-spin 
transformation by the matrix Urq. The transfer ma- 
trices are 'compressed' by the RG transformation 

(^RgA^Rg) /II^RgA^RgI^ — y 

(4g^rg)/I|4g^rgII 2 K q , (12) 

(^rg^s^rg) /II^rg^^rgI^ — > C q ■ 

We introduced the normalization factor || • || 2 in order 
to avoid the numerical over-flow in the expression of the 
partition function. 

The central issue concerns the definition of the RG 
transformation. In the density-matrix renormalization 
scheme, Ub.g is created by diagonalization of the reduced 
density matrix p which may be represented in a non- 
Hermitian (asymmetric) form 

p = Tr env \ip)(4>\. (13) 

The trace is taken over the spin variables bel ong ing to 
the environment as proposed by DMRG [l5|, l2fj | The 
states \tp) and \<j>) correspond to two parts of the whole 
lattice. The Boltzmann weight for these two parts can be 
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calculated as the product of the corner transfer matrices 
tf(<«s) = E CteatiWrtib) ■ • •C(^fc^), (14) 

4>(rtyts) =J2 C(^i)C(^iC 2 ) • • -C(K^) , (15) 

where we introduced the condition fc + ^ + 2 = g. The 
most optimal choice is to consider k = n q+ ^ and £ = ri g+ 4 
with n q given by Eq. (fTTj) . We have used letter /i for the 
2-state variable of 4> just for distinction from a of tp, and 
this choice is convenient when we construct the reduced 
density matrix. The normalized partition function can be 
written as -2(3, 9 ) = ("0 !</>)■ As an example, we obtain k = 
3 and I — 2 for the (3, 7) lattice with the corresponding 
Boltzmann weights 

Vteafc) = E cKaa)CKi6)CK 2 6)CK 3 ^) , 

Si 6 

Notice that if q is even, fc = £ = I — 1 resulting in 
\ip) = \4>)- For this choice the reduced density matrix p is 
always Hermitian (symmetric). However, for any odd q, p 
becomes non-Hermitian (asymmetric). This may lead to 
severe numerical instabilities. In order to avoid them, we 
symmetrize the reduced density matrix. We, therefore, 
consider an equally weighted reduced density matrix 

+ ^(^7)^(^7)- (17) 

Having tested both formulations of the reduced density 
matrix, Eqs. (1131) as well as (|17[) . we encountered numer- 
ical instabilities for the non-Hermitian case only, espe- 
cially in the vicinity of the phase transition. Otherwise, 
both density matrix formulations yield the identical ther- 
modynamic properties. 

IV. MAGNETIZATION AND ENERGY 

Since the detailed analysis of the phase transitions deep 
inside the hyperbolic lattice is of our interest, we concen- 
trate on the bulk properties of a sufficiently large inner 
region of the lattice [a, HH although the influence of the 
system boundary is not negligible at all for the discus- 
sion of the thermodynamic properties of the system. The 
bulk spontaneous magnetization is an example where the 
value can be calculated by 

M = Tr (up) /Tr p (18) 

in the CTMRG formulation. Without loss of generality, 
we set the coupling constant J and the Boltzmann con- 
stant &b to unity, and all thermodynamic functions are 
evaluated in the unit of fee- 
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2 4 6 8 10 12 14 16 18 20 



T 

FIG. 5: Spontaneous magnetizations M with respect to tem- 
perature T for 6 < q < 20. The inset shows the linear behav- 
ior of the cubic power of the induced magnetization M 3 with 
respect the magnetic field h around the transition tempera- 
turesT^. 



We now consider one-point functions of the Ising model 
on the series of (3,p) lattices in the thermodynamic limit. 
First of all, let us check the validity of our numerical pro- 
cedure as explained in the previous section. We perform 
a test calculation for the flat (3, 6) lattice. Keeping only 
m = 20 states of the multi-spin variables £ 0, O, [2(| , the 
obtained spontaneous magnetization is shown in Fig. [5] 
The estimated transition temperature T c = 3.641 is quite 
close to the exact value T c = 4/ In 3 « 3.64096 Q. 

Now, we focus on the hyperbolic surfaces. In Fig. [5] we 
also plot the temperature dependence of the spontaneous 
magnetization M for the coordination numbers from 7 < 
q < 20. The full and the dashed curves, respectively, 
distinguish the even and odd values of q. As we show 
later, the system is always off-critical whenever q > 7, 
even at the transition temperature. We, therefore, use 
the notation instead of T c (?) for q > 7; we also use 
Tpjp for q = 6 in order to unify the notation. 

If a small magnetic field h is applied at the transi- 
tion temperature Tp1~ 7 \ the cubed induced magnetiza- 
tion M 3 is always linear around h — 0. Thus, the model 
satisfies the scaling relation M(h,T pt ) oc h 1 / 8 with the 
scaling exponent 5 = 3. This value is known for the 
mean-field universality of the Ising model and is in full 
agreement with our previous results for the hyperbolic 
b > 5,4) lattices [TJ]- 

In order to observe the scaling relation of the spon- 
taneous magnetization M in a unified manner, we plot 
the squared magnetization M 2 in Fig. [5] with respect 

to the rescaled temperature by T^- 7 \ Near the point 
T = T { v l- 7) the mean-field behavior M(h = 0,T) oc 
{T^f - Tf with ft = i is detected. Note that on the 
(3, 6) lattice the exponent is (3 = | as displayed in the 
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FIG. 6: The squared spontaneous magnetization M 2 is linear 
with respect to the normalized temperature near the transi- 
tion point. This corresponds to = |. Inset: the linearity of 
the M 8 is observed only when q — 6 where f3 = ^. 




]n[T ( f-T] 



FIG. 7: Convergence rate of the effective scaling exponent. 
Inset: scaling of the phase transition temperatures T^p versus 
the integer q. 



inset. To detect the scaling exponent f3 in a more precise 
manner, we calculate the effective exponent 



Pee{T) 



d In 


M 


[h = 0,T<T^j 




d In 







(19) 



by means of the numerical derivative. The convergence 
of j3 e s(T) with respect to T$ - T is shown in Fig. It 
is apparent that the mean-field value j3 = i is detected 
for any q > 7, whereas we confirm /3 = ^ on the flat 
(3, 6) lattice only. The linear increase of the transition 




FIG. 8: Specific heat as a function of temperature. The open 
circles connected by the vertical dotted lines show the discon- 
tinuity. Inset: temperature dependence of the internal energy. 



temperature T^- 7 ^ with respect to q is shown in the inset 
where the linearity appears already around q > 8. This 
agrees with an intuition where the mean-field behavior 
becomes dominant for large coordination numbers. 

Let us analyze the specific heat (or the heat capacity) 
per bond 



dT 



(20) 



where E int is the internal energy per bond, or equiva- 
lently, the correlation function between the two nearest- 
neighbor spins 



-J((Ti(T i+ l) = -JTr((Ti<T l+1 p) 



(21) 



with Oi and <7i+i located at the center of the lattice. Fig- 
ure [5] shows the results for C v and E int . The internal 
energy E int is continuous for all the cases we computed. 
The presence of the kink in E int at the transition temper- 
ature for each q > 7 corresponds to the discontinuity in 
C v [H,[l^]. For these cases the scaling exponent a, which 

appears in the relation C v (h — 0, T) cx \Tpp —T\ , is 
zero. It is instructive to point out that both C v and E- lnt 
in the paramagnetic region are almost independent on q\ 
the tiny differences are hardly visible on the scale in the 
figure. 



V. ENTROPY AND CORRELATION 

Whenever the reduced density matrix p is defined, the 
von Neumann (or entanglement) entropy plj 



S=-Tr(plog 2 p) 



2^ log 2 Ui 



(22) 



can be used as a characteristic quantity which is of use for 
the classification of the phase transition. Figure [S] shows 
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FIG. 9: Temperature dependence of the von Neumann entan- 
glement entropy. The inset displays the dominant behavior 
of S for the (3, 6) lattice. 




FIG. 10: Decay of the density matrix spectra for the (3, 6) 
lattice (filled symbols) and the (3,9) lattice (open symbols). 



the temperature dependence of S which remains finite for 
q > 7 even at the transition temperature. The entropies 
in the paramagnetic region are also almost independent 
on q if q > 7 as found for C v and Ei nt . 

The decay rate of the density matrix eigenvalues uj i 
is shown in Fig. [10] on a semilogarithmic scale for both 
(3, 6) and (3, 9) lattices. We confirm a power-law decay 
in uj i only at the transition point of the (3, 6) lattice. 
Note that the eigenvalues uj i decrease exponentially for 
q > 7 at the transition temperature. 

The exponential decay of the density matrix spectra is 
also reflected in the correlation function 

G it j = Tr {ui<Jjp) (23) 

between two distant sites i and j. We place the spin <Ji at 




FIG. 11: Decay of the correlation functions with respect to the 
distance \i — j\. The open and the full symbols, respectively, 
correspond to the (3,6) lattice calculated at T = 3.0, 3.641, 
and 5.0, and the (3, 9) lattice at T = 6.0, 7.608, and 9.0. 



the center of the system and <jj at the system boundary. 
Therefore, as the lattice expands its size via the recursive 
steps in CTMRG, the distance between these two spins 
increases progressively. 

Figure [TT] depicts log 10 (Gi t j) as a function of \i — j\ 
for the (3, 6) lattice (open symbols) and the (3, 9) lattice 
(full symbols) . It is evident that the correlation functions 
always decay exponentially on the (3, 9) lattice regardless 
of the temperature. We remark that an analogous expo- 
nential decay of Gij has been observed for all q > 7 (not 
shown). On the (3,6) lattice, the correlation function 
decays as a power law at the transition temperature, as 
seen in the inset. 

In the following, we compare the Gaussian curvature 
associated to the (3, q) lattice with the correlation length 
at the transition temperature. There are several ways to 
define the correlation length £ q [l4|, H3|. For example, 
the decay rate of the correlation function directly pro- 
vides £ 9 . This is straightforward, but the region of the 
distance for the fitting analysis has to be valued care- 
fully. Another possibility consists in using the largest 
eigenvalue Ao (q) and the second largest one Ai (q) of the 
row-to-row transfer matrix where ^ q is determined from 
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In 



Ao(<?) 



Ai(g) 



(24) 



The relation can be generalized to the (3, q > 7) lattices, 
in analogy to our previous formulations for the (5, 4) lat- 
tice [24]], via the construction of the row-to-row transfer 
matrix 

T q (£ia a &\tiW a Q = £q(< 7 Wl£.lVa)£ q {cra£,2£,Wa) ■ ( 2 5) 

Using the notation of the recurrence scheme introduced 
in the previous section, we calculate £ 9 by use of Eq. (24) . 
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FIG. 12: Gaussian curvature K q with respect to the shifted FIG. 13: Asymptotic scaling of the correlation length £ q at 



phase transition temperatures for 6 < q < 20. The inset 



Eq. (I26|l while that on the right shows the correlation length 
in the vicinity of the phase transition. 



the transition temperatures T^' with respect to R q . The 
thin dotted lines serve as a guide to the eye. The inset shows 
the difference A^e-q(Tp^) in Eq. (|27[) with respect to q on a 
double logarithmic scale. 



The Gaussian curvature K q that corresponds to (3,q) 
lattice is given by [22j 



K q = 



(iR q y 



= — 4 arccosh 



2 sin 



(26) 



where R q is the curvature radius of the hyperbolic sur- 
face. Recall that K q must be zero on the Euclidean flat 
space (q = 6). Figure [121 shows the relation between K q 

and the shifted transition temperature T^f 1 — . The 
lower-left inset shows complementary information about 
R q . The correlation function ^ q calculated around the 
phase transition for three different g's is plotted in the 
upper- right inset. Notice that £ q reaches its maximum 
at the phase transition which is not well visible as q in- 
creases. 

Figure [T3] shows the dependence of the correlation 
length £ g (T) at the transition temperature with respect 
to the curvature radius R q . In order to collect these 
data, we performed extensive calculations up to 32 dig- 
its numerical precision for the value of q as large as 
q = 10 000 000 where the corresponding Gaussian cur- 
vature K W 7 is approximately 900. Note that both quan- 
tities diverge on the (3,6) lattice, and therefore ^(X^?) 
and Rq are not shown. Tet us focus on the limit R q — > 
which corresponds to q — ¥ oo. Evidently, the correlation 
length £ g decreases to zero as q tends toward infinity (the 
circles). Applying a least-square fit, we obtain the rela- 
tion £ q = lA4(iR q )°- 90S as shown by the thick dot-dashed 
curve. If we consider the error in the calculation of the 
correlation length, we can conjecture that £ g is propor- 
tional to R q . 

Recall that the specific heat C v , the internal energy 
£i n t, and the entanglement entropy 5* turned out to be 



weakly dependent on the value of q in the paramagnetic 
region T > T$ for q > 6. Thus, it can be conjectured 
that the disordered state is not modified by the presence 
of the negative curvature. We, therefore, compare £ 9 >7 
just at the transition temperature T$ with the corrc- 



T$. These 



lation length £ 6 at the temperatures T 

values are plotted in Fig. [13] by the asterisks. Since 
almost linearly increases with q for large values of q, the 
dotted line goes to the origin of the graph. The circles 
and the asterisks in Fig. [T3] are of the same order for 
all q, and this fact supports our conjecture that R q rep- 
resents the only characteristic length of the hyperbolic 
lattice and that the phase transition occurs at the tem- 
perature where £ q is of the same order as R q . Note that 

CeC^) > CgCT$?) is always fulfilled as plotted in the 
inset of Fig. [T3] where we show the difference 



A^-rf) 



e6«)-e,W) 



(27) 



The relation £, 6 (T^ ] ) > £q{T$) may be explained by the 
effect of the negative curvature that prevents from a kind 
of loop-back of the correlation effect. Such suppression is 
also expected to be present in higher-dimensional hyper- 
bolic lattices and could be analytically studied by means 
of the high temperature expansion. 

We conjecture the reason why the correlation length 
remains finite even at the phase transition temperature 
T^p for q > 6, as follows. First of all, the hyperbolic 
plane contains the typical length scale R q , and it might 
prevent scale invariance of the state expected at the crit- 
icality. A more constructive interpretation could be ob- 
tained from the observation on the row-to-row transfer 
matrix. The calculation of £ q by means of Eq. (PHI) re- 



8 



quires diagonalization of the row-to-row transfer matrix 
Tqi^i^a^Wi^'a^) m Eq. (j2"5)l . The matrix corresponds to 
an area which connects (transfers) the row of the neigh- 
boring spins {£iCTa<?2} with the adjacent ones {Cl^^l- 
The shape of this area is very different from the stan- 
dard transfer matrix on the Euclidean lattice, which cor- 
responds to a stripe of constant width. On the hyper- 
bolic surfaces, however, this distance between the spin 
rows is not uniform. The distance is minimal at the cen- 
ter of the transfer matrix, i.e., between the two spins 
a a and a' a , and it increases exponentially with respect 
to the deviation from the center to the direction of spin 
rows. Such a geometry [24{ could be imagined from the 
recurrence construction in Eq. ([§]). As a consequence, the 
transfer matrix has an effective width, which is of the or- 
der of the curvature radius R q . The region outside this 
width contributes as a sort of the boundary spins that im- 
poses mean- field effect to the bulk part. This situation is 
analogous to the Bethe lattice, being interpreted here as 
(oo, q)-lattices. [3]. Thus the Ising universality could be 
observed only when the correlation length £ g is far less 
than the curvature radius, f 5 <®C R q . As the length £ g 
increases toward the transition temperature, we expect 
a transient behavior to the mean-field behavior around 
the point when £ g becomes comparable to R q . We are 
confirming these conjectures and the details would be re- 
ported in our subsequent work. 

VI. CONCLUSIONS 

We have presented a detailed analysis of various non- 
Euclidean lattices forming surfaces with hyperbolic cur- 
vatures. In addition to our previous works on the (j>, 4) 
lattices, we studied the complementary situation repre- 
sented by the (3, q) lattices. This task required a reformu- 
lation of the existing CTMRG algorithm. We, therefore, 
considered the half-row transfer matrices and the corner 
transfer matrices including asymmetric (non-Hermitian) 



cases. For the lattices with odd q's, we symmetrized the 
density matrix by the way which has been accepted by 
the DMRG community H(|. 

We treated the Ising model on the (3, q) lattice with 
coordination number 6 < q < 10 7 . The phase transi- 
tion temperatures are determined from the analysis of 
the magnetization, internal energy, specific heat, and the 
von Neumann entanglement entropy. We have shown 
that the transition temperature T^ 1 linearly increases 
with q for larger values of q. The scaling behavior of the 
thermodynamic functions, including their related scaling 
exponents a — 0, /3 — ^, and 6 — 3, obeys the mean-field 
universality class. The mean-field nature of the hyper- 
bolic surfaces is characterized by the exponential decay 
of the reduced density matrix eigenvalues and the corre- 
lation functions even at the transition temperature. 

We further evaluated the radius of the Gaussian cur- 
vature R q for the generic (3, q > 6) lattice geometry and 
compare it to the results for the correlation length ex- 
tracted from the row-to-row transfer matrix. We found 
a strongly suppressed correlation length £ 9 < 1 at the 
transition point for any q > 7. We conjecture that ^ q is 
proportional to R q in the large q limit. 

In order to elucidate the origin of the mean-field uni- 
versality induced by the hyperbolic geometry, our future 
studies aim at the treatment of specific hyperbolic ge- 
ometries with non-constant Gaussian curvatures in order 
to systematically approach the Euclidean (flat) geometry. 
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